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Abstract 

We determine the asymptotic density 5k of the set of ordered fc-tuples (n±, ...,nfc) £ 
N fc , k > 2, such that there exists no prime power p a , a > 1, appearing in the canonical 
factorization of each m, 1 < i < k, and deduce asymptotic formulae with error terms 
regarding this problem and analogous ones. We give numerical approximations of the 
constants S k and improve the error term of formula (1.2) due to E. Cohen. 

We point out that our treatment, based on certain inversion functions, is applicable 
also in case k = 1 in order to establish asymptotic formulae with error terms regarding 
the densities of subsets of N with additional multiplicative properties. These generalize 
an often cited result of G. J. Rieger. 
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1 Introduction 

Let k > 2 be a fixed integer. What is the asymptotic density 5 k of the set of ordered 
A;-tuples (m, —,n k ) £ N fc , such that there exists no prime power p a , a > 1, appearing in 
the canonical factorization of each nj, 1 < i < k ? 

This problem is analogous to the following one: What is the asymptotic density of 
the set of A:-tuples which are relatively prime, i.e. /c-tuples (n\, ...,njfc) € N fe such that there 
exists no prime p, appearing in the canonical factorization of each m, 1 < i < k ? 

It is known that dk = 1/C(k), where ^ is the Riemann zeta function, and this value can 
be considered as the probability that k integers (k > 2) chosen at random are relatively 
prime. More precisely, 

(1.1) N k (x) := #{(m,...,n k ) £ (Nn [l,x}) k : gcd(m, ...,n k ) = 1} = -^x k + R k (x), 

where R k (x) = 0(x k ~ 1 ) for k > 2, R%(x) = 0{x\ogx) for k = 2, and d k = lim^.+oo N k (x)/x k 
= l/C(k). This result goes back to the work of J. J. Sylvester [9] and D. N. Lehmer 
[3], see also J. E. Nymann [5]. 
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There are several generalizations of (1.1) in the literature. For example, let S be an 
arbitrary subset of N. Then 



(1.2) N k (x,S) := n fc ) € (Nn[l,i]f : gcd(n u ...,n k ) £ S} = ^-l x k + T k (x), 



and Tfc(ar) = 0(x fc x ) for k > 2, T2(x) = 0(xlog 2 x) for k = 2, for every 5CN, due to E. 
Cohen [1]. Therefore the asymptotic density of the set of ordered fc-tuples (ni,...,n k ) for 
which gcd(ni, ■ ■■,n k ) belongs to 5 is lim^oo 7V fc (x, S)/x k = 

J. E. Nymann [6] shows that if the characteristic function of 7^ S 1 C N is completely 
multiplicative and if #{n : n £ S l~l [1, x]} = Ax + O(l), where A is the asymptotic density 
of S, then 



(1.3) #{(m, n fc ) G(Sf1 [1, x]) fc : 5 cd(m, n fc ) = 1} = A fc JJ 1 - - ) x k + R k (x), 



where R k {x) is the same as above. Therefore, if P k {n) denotes the probability that k 
integers (k > 2) chosen at random from S n [1, n] are relatively prime, then 



This result can be applied for S = {n : gcd(n,p\ ■ ■ -p r ) = 1}, where {pi, ...,p r } is a given 
finite set of distinct primes. 

Now return to the problem at the beginning. It is obvious that 5 k > d k = 1/Q{k) for 
every k > 2 and thus lim^oo 5 k = 1. Which is the exact value of 5 k ? 

In order to solve this problem we use the concept of the unitary divisor. For d, n 6 N, 
d is called a unitary divisor (or block divisor) of n if d\n and gcd(d,n/d) = 1, notation 
d\\n. Various other problems concerning unitary divisors, including properties of arithmeti- 
cal functions and arithmetical convolutions defined by unitary divisors, have been studied 
extensively in the literature, see for example [4] and its bibliography. Denote the greatest 
common unitary divisor of ni, ...,n k by gcud(n\, ...,n k ). 

Our question can be reformulated in this way: What is the asymptotic density 5 k of 
the set of ordered fe-tuples (ni,...,n k ) such that gcud(ni, ...,n k ) = 1, or more generally, 
gcud(ni, ...,n k ) £ 5 ? 

Furthermore, what is the probability that for k integers n\, n k chosen at random from 
S fl [1, n] one has gcud(rii, n k ) = 1 ? 

In this paper we determine the value 5 k and deduce asymptotic formulae with error ter- 
mes analogous to (1.1) -(1.3), regarding these problems. We give numerical approximations 
of the constants 5 k and also improve the error term of (1.2) of E. Cohen. 

The treatment we use is based on on the inversion functions jjl* s and us attached to 
the subset S. We point out that this is applicable also in case k = 1 in order to establish 
asymptotics regarding the densities of certain subsets S of N, generalizing in this way an 
often cited result of G. J. Rieger [7]. 



where 



00 

n=l 
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Note that the value <5 2 is given by D. Suryanarayana and M. V. Subbarao [8], 
Corollary 3.6.3, applying other arguments as those of the present paper. 

Using the concept of regular cross-convolution, see [11], [12], it is possible to deduce 
more general results, including (1.1) - (1.3) and (2.1) and (2.4) of this paper. We do not go 
into details. 



2 Results 

Let S C N. We say that S is (completely) multiplicative if 1 G S and its characteristic 
function ps(n) is (completely) multiplicative. Define the function p* s (n) by 

X>s(d) = Ps(n), neN, 

d\\n 

that is 

/fsrfa) = S Ps(d)v*(n/d), n G N, 

d||n 

where the sums are extended over the unitary divisors of n and //*(n) := //^(n) = (— l) w ( n ), 
w(n) denoting the number of distinct prime factors of n. 

Furthermore, let <j)[n) and 6{n) denote Euler's function and the number of squarefree 
divisors of n, respectively. 

Theorem 2.1 If k > 2 and S is an arbitrary subset o/N, then 

(2.1) #{(m, n fc ) G (N n [1, x]) fe : gcud( ni ,..., n k ) £ 5} = <5 fe (S)x fe + Vfc(x, 5), 

w/iere 

4(5) 

n=l 

and t/ie remainder term can be evaluated as follows: 

(1) Vk(x, S) = 0(x k ~ 1 ) for k > 2 and /or an arbitrary S, 

(2) V2(x,S*) = 0(xlog x) for an arbitrary S, 

(3) V2(x,S) = 0(xlog 2 x) for an S such that J2 n eS < 00 particular for every 
finite S ) and for every multiplicative S, 

(4) 1/2(2;, S) = 0(x) for every multiplicative S such that J2 P ^s | < 00 particular if 
the set {p : p ^ S} is finite). 

If S is multiplicative, then 

s k (s) = n 

v 

IfS = {l}, then 
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Theorem 2.2 If k > 2 and S is an arbitrary subset o/N, then the asymptotic densities of 
the sets of ordered k-tuples (n±, n^) swc/i i/tai gcud(n\, n k ) G S 1 and gcud(n±, n k ) = 
1 are Ofc(S') ana 1 6fc, respectively, given in Theorem 2.1. 



Theorem 2.3 Z-ei p n denote the n-th prime and for r e N let N = 10 r /2. T/ien 




(Pn ~ 1)" 



) 



n=l \ 



is an approximation of S k with r exact decimals. 

In particular, S 2 ~ 0.8073, 5 3 « 0.9637, <5 4 « 0.9924, o 5 « 0.9983, 5 6 « 0.9996, <5 7 « 
0.9999, wii/i r = 4 exarf decimals. 

Theorem 2.4 For k = 2 the error term B,2(x) of (1.2) can 6e improved into R(x, S), where 

(i) i?(x, 5) = O(xlogx) /or an S such that J2neS \ < 00 particular for every finite 
S ) and for every multiplicative S, 

(ii) R(x,S) = 0(x) for every multiplicative S such that J2 P <fs \ < 00 ^ n particular if 
the set {p : p ^ S} is finite). 

Remark. It is noted in [1] that if k = 2 and if the function fis(n) is bounded, cf. proof 
of Theorem 2.4 of the present paper, then the error term is R.2(x) = O(xlogx). 

Theorem 2.5 Suppose that S C N is multiplicative and min{a : p a £ S} > r > 2 for every 
prime p. Then 



is the asymptotic density of S. 

Remark. In the special case S = the set of If -void integers we reobtain from (2.2) 
the result of G. J. Rieger [7]. The if- void integers are defined as follows. Let if be a 
nonempty subset of N \ {1}. The number n is called if -void if n = 1 or n > 1 and there is 
no prime power p a , with a G if , appearing in the canonical factorization of n. 

Does the density exist for an arbitrary multiplicative subset S ? Yes, and it is d(S) 
given by (2.3), where the infinite product is considered to be when it diverges (if and 
only if J2 p <£s p = °°)- This follows from a well-known result of E. WlRSlNG [13] concerning 
the mean- values of certain multiplicative functions /. A short direct proof for the case / 
multiplicative and < f(n) < 1 for n > 1, hence applicable for the characteristic function 
of an arbitarary multiplicative S, is given in the book of G. Tenenbaum, [10], p. 48. 

Theorem 2.6 let k > 2 and suppose that S is a completely multiplicative subset o/N such 



that #{n : n G S n [1, x]} = Ax + 0(1). Then 

(2.4) #{(m, ...,n k ) e (Sn[l,x]) k : gcud( ni ,...,n k ) = 1} = A k p k (S)x k + T k (x), 



(2.2) 



d(S)x + 0«/x~), 



n<x 



where 



(2.3) 
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where 



and T fc (x) = 0{x k ~ 1 ) for k > 2, T 2 (x) = 0{x\og z x) for k = 2. 

If Qf(n) denotes the probability that for k integers n\,...,n k chosen at random from 
Sn[l,n] one has gcud(n\, ...,n k ) = 1, then 



lim Q$(n) = p k (S). 

n— >oo 

3 Proofs 

Proof of Theorem 2.1 Using the definition of /j,* s , the fact that d\\gcud(ni, ...,n k ) if and 
only if d||nj for every 1 < i < k, which can be checked easily, and the well-known estimate 



gcd( i n,m)=l 



which holds uniformly for x > 1 and m £ N, we obtain 

#{(ni,...,n fc ) G (Nn : gcud(n l7 ...,n k ) 6 5}= E p s {gcud{n u ...,n k )) = 

m,...,n k <x 

= E E E E A*S(d) = 

»H,...,n fc <xd||(ni,...,ra fc ) ni,...,n fc <xd||ni,...,d||n fc 

= E/4W E i = E^) E 1 

cKx aj<a:/d d<X \ a<x / d 

( 0i ,d) = l \(a,d) = l / 

l<i<fc 



d>x 



The given error term yields now from the next statements: 
(a) 

0(n) = (0{log 2 x), 8 = 1, 
£i n s 10(1), a >l. 



n<x 

32 



v fl 2 (n) _ fO(log 4 x), 8 = 1, 
^ n* -\0(1), a >l, 



El , 1 s 6»(n) _ logx ! 



„n x° '■ rrl n £ 
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(b) For an arbitrary S C N and for every n £ N, |//g(n)| < J2d\\nPs(d) < 
|/4(n)|0(n) < E d \\ n Ps(d)e(d)6(n/d) and 

v |M|Hjg(n) < ps(d)6(d) T 0(e) _ 

n<x d<£ e<x/d 

= (log 2 x V PsWM ) = { ^(log 2 x), if E~ x ^« < oo, 
I d / I Oflo^x). otherwise. 



d<x 



(c) If S 1 is multiplicative, then [i* s is multiplicative too, fi* s (p a ) = ps{p a ) — 1 for every 
prime power p a (a > 1) and fJ-*s(n) £ { — 1,0, 1} for each n £ N. 

(d) Suppose S is multiplicative. Then 

\Ps(p a Mp a ) _ 2 y y i-ps(p a ) < 

P a=l ^ p a=l ^ 

oo 



< 4 IE J 2 + E-)< o ° if E-< 



oo. 



It follows that in this case the series J2™=i is convergent. 

If S is multiplicative, then the series giving Sk(S) can be expanded into an infinite 
product of Euler-type. 

Proof of Theorem 2.2 This is a direct consequence of Theorem 2.1. 

Proof of Theorem 2.3 Consider the series of positive terms 

(p-i) 1 * 1 



where p n denotes the n-th prime. 

The iV-th order error of this series can be evaluated as follows: 

°° ( (Pn~ l) fe \ °o (p n — l) fe 

^ := „=¥ + i log V 1 + j*(pfe-i)-(P»-i) fc J < J£ +1 p£(p£ - 1) - (p» - !) fc < 

oo ^ oo ^ 

< E pfc _ 1 - E „2 _ X ' 

Now using that p n > 2n, valid for n > 5, we have 

W< V^™ 2 " 1 2 4li\2n-l 2n + l7 2(2iV + l)' 
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In order to obtain an approximation with r exact decimals we use the condition 

1 1 



2(27V + 1) - 2 



C - ■ 10 



-r 



and obtain N > ±(10 r - 1). 

The numerical values were obtained using the software package MAPLE. 
Proof of Theorem 2.4 Define the function ps{ n ) by 

^Hs{d) = ps{n), n G N, 

d\n 

that is 

Ps(n) = Y Ps{d)n{n/d), n G N, 

d\n 

where p{n) := ^{^(n) is the Mobius function, see [1]. We have 

N k (x,S) := #{(m,...,n fc ) G (Nn[l,x]) fc : <7cd(m, n fc ) 6 5} = 

= E p s (gcd(m,...,n k ))= Y Y ^( d ) 

m,...,n k <x m,...,n k <x d\(m,...,n k ) 

and the proof runs parallel to the proof of Theorem 2.1. 
Proof of Theorem 2.5 

tfi(s, 5) = ^ Ps(n) = EE M<0 = * E ^ + °(E 

n<x n<x d\n d<x d<x 

Here p$ ls multiplicative, Ps(p a ) = Ps(p a ) ~ Ps(p a l )> a ^ 1 an d since p,p 2 , ...,p r ~ 1 G S 1 
we have ps{p) = Ps{p 2 ) = ••• = Psip 7 "^ 1 ) = for every prime p. Hence for each n G N, 
l/^sWI < PL r { n )i where L r is the set of r-full numbers, i. e. L r = {l}U{n > 1 : p\n p r \n}. 
We get 

mix, s) = d(s)x + o(x Y Eh j^) + o(E (<*)), 

and using the elementary estimate 

^ /9Lr (n) = C v ^ + 0(^), 

n<a; 

where C is a positive constant, due to P. Erdos and G. Szeker.es [2], obtain the given 
result. 

Proof of Theorem 2.6 

.,n k ) G (Sn [l,x]) k : gcud(ni,...,n k ) = 1} = 
= E Ps( n i)- E ^s(njfc) Y = E P«( n i)- E Ps( n k) E = 

ni<x n k <x d\\(n\,...,n k ) n\<x n k <x d||ni,...,d| |n fc 



E^*( d ) E Ps(dai)... Y Ps(da k ) = 

d<X ai<x/d a k <x/d 

(oi,d) = l (a k ,d) = l 



7 



d<x 



a<x/d 
\ (a,d)=l 



J 



Here we use the estimate, valid for every leN, 



Ps{n) = Y Ps{n) P( d ) = 

n<x d\gcd(n,£) 



gcd(n,£) = l 



d\l 



e<x/d 



d\t 



= AxH(l-±) + 0(e(£)) 



pes 



and obtain the desired result, see the proof of Theorem 2.1. 
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